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Abstract 

In the framework of the semihard (kx factorization) QCD approach, we consider the pho- 
toproduction of D*^ mesons associated with two hadron jets and the D*^ production in 
DIS at THERA conditions with the emphasis on the BFKL and CCFM dynamics of gluon 
distributions. In the photoproduction of D*^ mesons the attention is focused on the variable 
z-y, which is the fraction of the photon momentum contributed to a pair of jets with largest 
Pt- We show that our theoretical results are sensitive to the BFKL type dynamics which 
may be investigated at THERA energies. We also discuss possible effect of J/tp meson spin 
alignement, which is thought to be a vivid manifestation of gluon off-shellness. 
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1. Introduction 



The experimental results on heavy flavour production processes obtained by the HI 
and ZEUS ||2|, ^ collaborations at HERA provide a strong impetus for further theoretical and 
experimental studies in a new energy region at THERA conditions. 

In due time, the experimental data have been compared with next-to-leading order (NLO) 
perturbative QCD calculations using the 'massive' and 'massless' schemes. The measured 
cross sections generally lie above the predicted level, and an agreement between the theoretical 
and experimental results can only be achieved using some extreme parameter values. In 
particular, the production rates of D*^ mesons in the NLO massive scheme @ require as 
low quark mass as rric = 1.2 GeV and as sharp charm fragmentation function as e = 0.02 
(in the Peterson parametrization) . However, even within this set of parameters, the shapes 
of the D*^ transverse momentum and rapidity distributions cannot be said well reproduced. 
A good agreement between the massless scheme |^] and the measured pt{D*) (though not 
r]{D*)) spectrum was achieved upon introducing an additional charm excitation contribution 
assuming an incredibly large charm content in the photon structure functions [Q: c{x) ~ u{x). 

The so called kx factorization, or the semihard approach (SHA) on one hand may 

give a reasonable solution for some of the above problems. On the other hand the significance 
of kT factorization (semihard) approach becomes more and more commonly recognized. Its 
applications to a variety of photo-, lepto- and hadroproduction processes are widely discussed 
in the literature |11| - |T^. In many cases a remarkable agreement is found between the 
data and the theoretical calculations regarding the photo- [16| and electroproduction |17| of 
D*^ mesons, forward jets |l^], as well as for specific kinematical correlations observed in 
the associated D*^ +jets photoproduction [19| at HERA and also the hadroproduction of 
beauty 21], Xc |22] and J/ip [23| at Tevatron. The theoretical predictions made in ref. 



|24| has triggered a dedicated experimental analysis |25] of the J/^l^ polarization (i.e., spin 
alignement) at HERA energies. 

To some extent, the SHA based on the BFKL [^] gluon dynamics includes the relevant 
effects of higher order contributions |27, 28|. It has been also demonstrated in [^] that 
the SHA effectively imitates the anomalous coupling of the resolved photon and an ad hoc 
contribution from the resolved photon is no longer needed at least for the description of the 
Xy distribution in D*^ photoproduction at HERA energies. 

In the present paper we use the semihard QCD approach to predict some features of the 
D*^ and J/ip production processes in the new energy region of THERA collider. 



2. The semihard QCD approach 



The production of J/ip mesons and open-flavoured cc pairs is described in terms of 
the photon-gluon fusion mechanism. A generalization of the usual parton model to the fcy- 
factorization approach implies two essential steps. These are the introduction of unintegrated 
gluon distributions and the modification of the gluon spin density matrix in the parton-level 
matrix elements. 

At first we consider the relevant partonic subprocesses. Let ki, k^, k^ and be the 
momenta of the initial state photon, the initial state gluon, the final state gluon and the 
final state J/ip, respectively, ei, €2, £3 and the polarization vectors, and k=ki+k2- The 
photon-gluon fusion matrix elements then read: 

Mi'yg tpg) = tr{4i {pc- fi\ + rric) ^ (- ^g- /;3 + mc) ^ J(e^)} 

x[A;i - 2(pcA:i)]~n^3 + 2(^5^:3)]"-^ + 5 permutations (1) 

Similarly, for the production of an open-flavoured cc pair (see fig. lb): 

A^(7g cc) = u{pc){/^i {pc- h + ^c) h [kf - 2A;iPc]~^ 

+ h {Pc- h + mc) h [ki - 2k2Pcr^}uipc) (2) 
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In the expression (||), the projection operator ||2^ J{^ip) = {/i^c + fric)/ guarantees 
the proper spin structure of the cc state, and the charmed quarks are assumed to each carry 
one half of the J/ip momentum, pc=Pc=P^/'^, W'c=m^/2. The formation of a meson from 
the cc pair is a nonperturbative process. Within the nonrelativistic approximation we are 
using, this probabiUty reduces to a single parameter related to the meson wave function at 
the origin |^(0)p, which is known for J /ijj and T families from the measured leptonic decay 
widths. 

The evaluation of traces in ([l|)-(^) is straightforward and is doneusing the algebraic 
manipulation system FORM ||3^. The complete set of matrix elements has been tested for 
gauge invariance by substituting the gluons momenta for their polarization vectors. 

The multiparticle phase space H d^Pi/i'^Ei) 5'^{J2Pin — J2Pout) is parametrized in terms 
of transverse momenta, rapidities and azimuthal angles: ^^^=^dpfrpdyi^ . Let (pi and 4>2 be 
the azimuthal angles of the initial photon and gluon, (p^, (p^, (pc and (pc the azimuthal angles 
of the partonic subprocess products (i.e., oi & J/ip and a gluon, or a charmed quark and 
an antiquark, respectively) and 1/3, y^, and T/g their rapidities. Then, the fully differential 
cross sections read: 

spins colours 



xTg{x2, A;2T, ^ ) dkiT dk2T dp^T dys dy^ — (3) 



27r 2tt 2tt 

\^ i \^ 

16-K X2 4 

d(pc 



daiep ^ c-cX) = 7 E ^ E - cc)\^ 



xj^g{x2, kjr^, fi^) dkfj. dklrp dplrp dye (4) 
The phase space physical boundary is determined by the inequality 

G{sXklklklkl)<Q, (5) 

where fcs and k^ denote the partonic subprocess final state momenta, s = {ki + ^2)^, t = 
{ki — k'if', and G is the standard kinematical function |31|. 

When calculating the spin average of the matrix element squared, we substitute the full 
lepton tensor for the photon polarization matrix: 



e'ter = ^namp^, - A{p,ki)g^^^]/ {k^Y (6) 

(including also the photon propagator factor and photon-lepton coupling). For the off-shell 
incoming gluon we take Q 

l^ = ki^Tk-2T/\k2T\^- (7) 

This formula converges to the usual ^ e^e*^ = —g^^ when kT — > 0. The final state gluon in (|l|) 
is assumed on-shell, ^ £363^ = —g^'^- The J/'ip polarization vector is defined as an explicit 
four- vector. In the frame with z-axis along the J/ip momentum, p^ = (0, 0, E^), it 
reads for different helicity states: 

eV('i=±i) = (1' 0' 0)/^> ^i>{h=o) = (0, 0, E^, \p^\)/m^. (8) 

The initial photon and gluon momentum fractions xi and X2 are calculated from the 
energy-momentum conservation laws in the light cone projections: 

{h + k2)E+pn = xi^/s = m3Texp(y3) + |m4T|exp(?/4), 

(9) 

{h + k2)E-pn = X2\fs = msT exp(-y3) + |m4r| exp(-y4). 
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The multidimensional integration in (Sh, (^) has been performed by means of Monte-Carlo 
technique, using the routine VEGAS | 3^ . 

Another important ingredient of the semihard approach is the so called unintegrated 
gluon distribution J^{x, kj^, fjp), which determines the probability to find a gluon carrying 
the longitudinal momentum fraction x and transverse momentum fcy- In calculations we use 
two different sets of unintegrated gluon distributions. One of them is based on the approach 
of ref.|Q and is costructed as a leading-order perturbative solution of the BFKL equations. 
Technically, the unintegrated gluon density J-g{x, k'^, fj?') is calculated as a convolution of the 
ordinary gluon density G{x,fj?) with universal weight factors: 



J' g [x, krp, n 



g{r],kT,n'^ 



T] rj 



(10) 



g{r],kT,^i^) 



■qkj, 

as 
7] k?p 



70(2^^0, ln(l/r?)ln(/.2/4)), 
10(2^/0, ln(l/r?)ln(4/^2))^ 



k^ < A*^, 



k^ > /i^. 



fill 



(12) 



where Jq and Iq stand for Bessel functions (of real and imaginary arguments, respectively), 
and Os = as/Sn. The LO GRV set |34] was used here for the boundary conditions. Another 
set of unintegrated gluon densities was extracted from a numerical simulution of the CCFM 
equations |35| and then tabulated in the form of a FORTRAN code [18]. Finally, the charm 
quarks were converted into D*^ mesons using the Peterson fragmentation function |3G]. 

In the paper we used the standard GRV parametrization |34] for the collinear gluon 
density, from which the unintegrated gluon distribution was developed according to eqs. (|lO|)- 
([T^). Some other essential parameters were chosen as follows: the charm quark mass nic = 
1.5 GeV, the Peterson fragmentation parameter e = 0.06, the overall c D* fragmentation 
probability 0.26. The Pomeron intercept A was regarded as free parameter, and then the 
value A = 0.35 has been extracted from a fit to the experimental pt{D*) spectrum measured 
by the ZEUS collaboration Q. Close estimates for A have also been obtained by many other 
authors, see, e.g. p7|, B^]. Since the agreement with the data achieved within this set of 



parameters was really good 17|, we continue using it in the present calculations. 



0,1 Q,? 0,3 0.^ O.f] 0,5 0.7 0.8 0.9 
Xgam 



0.1 0.2 0.0 0.^ 0..="] 0.6 0.7 0.3 0.9 



Figure 1: The differential cross section da/dx^ for Q'^ < I GeV^ with BFKL and CCFM 
unintegrated gluon distributions at THERA. 
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3. Numerical results 



3.1 D*^ and dijet associated photoproduction at THERA 

The ZEUS collaboration has measured the associated charm and dijet production [|| as 
a further test of the underlying parton dynamics. In these measurements, the quantity of 
interest is the fraction of the photon momentum contributing to the production of two jets 
with highest Et, which is experimentally defined as 

= [EiT exp(-r/i) + E2T exp(-r/2)]/(2Se y) (13) 
with EiT and 7]i being the transverse energy and rapidity of these hardest jets. 



In the ref. |17] the theoretical calculations have been made within the SHA with differ- 
ent unintegrated gluon distributions at HERA energies. In Fig. 1 we present the results of 
similar theoretical calculations made within the semihard approach with BFKL and CCFM 
unintegrated gluon distributions at THERA energies. The simulation procedure consists in 
generating a photon-gluon fusion event using the off-shell matrix elements and the uninte- 
grated gluon distribution functions described in Section 2. 

The basic 2^2 partonic subprocess gives rise to two high-energy quarks, which can 
further evolve into hadron jets. Actually, as the matter of some reasonable approximation, 
the calculations were restricted to parton level, and so the produced quarks (with their known 
kinematical parameters) were taken to play the role of the final jets: ET{jeti^2) = ET{q,q)- 

The two quarks are accompanied by a number of gluons radiated in the course of the 
gluon evolution. It has been mentioned already that, on the average, the gluon transverse 
momentum decreases from the hard interaction block towards the proton. As an approxima- 
tion, we assume that the gluon closest to the quark block with its momentum k' compensates 
the whole transverse momentum of the virtual gluon participating in the hard interaction: 
k'x — —kT, while all the other emitted gluons are collected together in the 'proton remnant', 
which is assumed to carry only a negligible transverse momentum compared to /c't- This 
gluon gives rise to a third hadron jet with Ej- = |A;'t|- 

From among the three hadron jets represented by the quark, antiquark and gluon we 
choose the two carrying the largest transverse energies, and then calculate the quantity 
according to its definition given by equ. ( p^ |^. 

In a significant fraction of events, the gluon radiated from the BFKL cascade appears to 



be harder than one or even both of the quarks produced in hard parton interaction |19]. In 
fact, the specified events are responsible for the wide plateau seen in the x^ distribution in 
Fig. 1. 

3.2. Deep inelastic D*^ production at THERA 

The process of deep inelastic D*^ production at HERA is truly semihard because of the 
presence of two large scales: the virtuality of the exchanged photon (Q^) and the charm mass 
(m^), both being much larger than Aqcd but much smaller than s. Therefore, experimental 
data concerning the D*^ production in DIS at THERA provide a strong impetus for further 
theoretical studies of this process. 

In Fig. 2 the theoretical predictions on the differential cross sections of deep inelastic 
D*^ production are shown for the THERA kinematical region: 1 < < 1000 GeV^, 
1.5 < pt{D*^) < 15 GeV and \ri{D*^)\ < 1.5. Different curves in Fig. 2 correspond to 
the BFKL and CCFM unintegrated gluon distributions. At HERA energies the SHA calcu- 



lations with BFKL unintegrated gluon distribution have shown |17] some shift to negative 
r]{D*) with respect to the ZEUS data. When we have used the CCFM unintegrated gluon 
density from MC generator CASCADE M] with JETSET based fragmentation function [p9|] 



^The full FORTRAN code is available from the authors on request. 
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Figure 2: Differential cross sections for deep inelastic D production with BFKL and 
CCFM unintegrated gluon distributions at THERA as functions of: (a) log^oQ^; (^) logio 
(c) W, (d) pt{D*), (e) r]{D*) and (/) z{D*). 

implemented in, we obtain good agreement between our theoretical results and the ZEUS 
experimental data @ also for da / dr]{D*) f[% . 

3.3. J/tp photoproduction at THERA 

The role of the gluon virtuality may be seen in Fig. 3, where we show the results of 
calculations for J/^l' photoproduction made with the BFKL unintegrated gluon distribution. 
The results correspond to THERA conditions, i.e. electron proton collisions at the energy 
^/s = 1000 GeV, where no other cuts were applied except the photoproduction limit < 1 
GeV^ and the inelasticity requirement 0.4 < z < 0.9. 

The effects of initial gluon off-shellness may be, best of all, seen in the transverse momen- 
tum spectra ||2^, because the gluon virtuality is proportional to its transverse momentum: 

= —kj,/{l—x). In contrast with the conventional (massless) parton model, the SHA shows 
that the fraction of J/ tp mesons in the helicity zero state increases with their transverse mo- 
mentum pt- a deviation from the parton model behaviour becomes well pronounced already 
from Pt > 3 GeV at HERA energies p^ , and at px > G GeV the helicity zero polarization 
tends to be even dominant (Fig. 3c). 

Qualitatively, the difference between the model predictions referes to different origins of 
the J/ip transverse momentum. In the case of conventional parton model J ftp meson obtains 
its transverse momentum from the hard photon gluon interaction, while in the SHA there is 
also a large the contribution from the initial gluon transverse momentum. 

The degree of spin alignement can be measured experimentally since the different polar- 
ization states of J/^p result in significantly different angular distributions of the J/ip ^ l~^l~ 
decay leptons: 

dTh=±i/dcose = l-cos'^@, dTh=o/dcose = l + cos^e (14) 
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Figure 3: Differential cross sections for inelastic J/^ photoproduction with BFKL uninte- 
grated gluon distributions at THERA: (a) Inclusive J/^ transverse momentum distribution, 
(b) the same, but for J/^ zero helicity states only, (c) the fraction of J/^ mesons in helicity 
zero state (degree of spin alignement). 

Here is the angle between the lepton and J/V"iirections, measured in the J/V' meson rest 
frame. Evidently the most informative regions relate to cosQ = ±1. 

4. Conclusions 

In the framework of semihard QCD approach, we obtained some predictions for the cross 
sections of inclusive D*^ meson production at THERA conditions using different unintegrated 
gluon distributions driven by the BFKL and CCFM evolution equations. 

We have considered the photoproduction of D*^ mesons associated with two hadron jets 
and also D*^ production in DIS at THERA conditions, which may be a sensitive indicator 
of the underlying parton dynamics. The results of the simulations show that theoretical 
results are very sensitive to BFKL type dynamics, in particular, to the unintegrated gluon 
distribution in the proton. 

We have considered also the effects of initial gluon off-shellness in SHA for J /ip meson 
photoproduction at THERA energies. Gluon virtuality connected with its transverse mo- 
mentum is one of the inherent properties of noncollinear (BFKL) parton evolution theory. 
Compared to traditional (collinear) parton model, gluons are characterized by a different 
spin density matrix. The latter is found to affect the polarization of J /ip mesons produced 
in ep collisions via photon gluon fusion subprocess. The effect is best pronounced at large 
J/il) transverse momenta and can be detected experimentally by measuring the J/if) ^ 
decay lepton angular distributions. We recommend the above process as a direct probe of 
the gluon virtuality, which can eventually testify for the validity of BFKL gluon evolution. 

Thus the experimental and theoretical investigations in the new kinematic region of 
THERA collider will provide additional tests of the semihard (fc-r factorization) approach 
and, in particular, of the "universality" of unintegrated gluon distribution. 
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